Influence of surface roughness on superhydrophobicity 
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Superhydrophobic surfaces, with liquid contact angle 6 greater than 150°, have important practical 
applications ranging from self-cleaning window glasses, paints, and fabrics to low-friction surfaces. 
Many biological surfaces, such as the lotus leaf, have hierarchically structured surface roughness 
which is optimized for superhydrophobicity through natural selection. Here we present a molecular 
dynamics study of liquid droplets in contact with self-afhne fractal surfaces. Our results indicate 
that the contact angle for nanodroplets depends strongly on the root-mean-square surface roughness 
amplitude but is nearly independent of the fractal dimension Df of the surface. 



The fascinating water repellents of many biological sur- 
faces, in particular plant leaves, have recently attracted 
great interest for fundamental research as well as prac- 
tical applications |1, _2, _3, 4_, _5, 6_, J^ &}■ The ability of 
these surfaces to make water beads off completely and 
thereby wash off contamination very effectively has been 
termed the Lotus effect, although it is observed not only 
on the leaves of Lotus plant, but also on many other 
plants such as strawberry, raspberry and so on. Water 
repellents are very important in many industrial and bi- 
ological processes, such as prevention of the adhesion of 
snow, rain drops and fog to antennas, self-cleaning win- 
dows and traffic indicators, low-friction surfaces and cell 
mobility HEIHI- 

Most leaves that exhibit strong hydrophobicity have hi- 
erarchical surface roughness with micro- and nanostruc- 
tures made of unwettable wax crystals. The roughness 
enhances the hydrophobic behavior, so that the water 
droplets on top tend to become nearly spherical. As a 
result the leaves have also a self-cleaning property: the 
rain drops roll away removing the contamination parti- 
cles from the surface, thanks to the small adhesion energy 
and the small contact area between the contaminant and 
the rough leafflj|. 

The hydrophobicity of solid surfaces is determined 
by both the chemical composition and the geometrical 
micro- or nanostructure of the surface |l2lll3lll4l |. Under- 
standing the wetting of corrugated and porous surfaces is 
a problem of long-standing interest in areas ranging from 
textile science jl5| to catalytic reaction engineering [ig. 
Renewed interest in this problem has been generated by 
the discoveries of surfaces with small scale corrugations 
that exhibit very large contact angles for water and other 
liquids — in some cases the contact angle is close to 180°. 
Such surfaces are referred to as superhydrophobic [l2J. 

The contact angle 9 between a flat solid surface and 
a liquid droplet depends on the relation between the 
interfacial free energies per unit area: solid/liquid 7si, 
solid/vapor 7sv and liquid/vapor 7iv. The Young's equa- 
tion 7si -|- 7ivCOS0 = 7sv, results from the minimization 
of the total free energy of the system on a flat substrate 
surface. Complete wetting corresponds to 6* = 0°, and 



typically happens on solids with high surface energy 7sv- 
Liquids on low energy surfaces tend to form droplets with 
high contact angle 9. 

It is well known that the roughness of a hydrophobic 
solid (with 9 > 90° on the flat substrate) enhances its 
hydrophobicity. If the contact angle of water on such 
flat solids is of the order of 100° to 120°, on a rough 
or microtextured surface it may be as high as 150° to 
175° |llLll8lll9l |. Two distinct models have been proposed 
to explain this effect. The Wenzel model 20] considers 
the increase of contact area due to the surface rough- 
ness: this leads to an increase of the effective free energy 
of the solid/liquid interface, making the surface more hy- 
drophobic. The contact angle 9o on the rough surface is 
obtained from the contact angle 9 on the microscopically 
flat surface of the same material through this equation 



cos9o = rcos9 (Wenzel model). 



(1) 



where r — A/Aq is the ratio between the real substrate 
area and the nominal (or projected) area Aq. 

The Cassie model _2l| assumes that some air remains 
trapped between the drop and the cavities of the rough 
surface. In this case the interface free energy 7si must 
be replaced by a weighted average of three interface free 
energies 7si, 7iv and 7sv, with the weights depending on 
the fraction (p of the area where the contact between the 
liquid and the solid happens. The contact angle is given 
by 



cos 6*0 



-1 + 0(1+' 



(Cassie model). (2) 



Quere states that there exists a critical contact angle 9c 
such that the Cassie state is favored when 9 is larger than 
9c ,.22j|. For a micro- or nano structured substrate, usu- 
ally the droplet stays in the Cassie state, but the Cassie 
state can switch (irreversibly) to the Wenzel state when 
the droplet is pressed against the substrate [23. The 
Wenzel droplets are highly pinned, and the transition 
from the Cassie to the Wenzel state results in the loss 
of the anti-adhesive properties generally associated with 
superhydrophobicity. 



Many surfaces in nature, e.g., surfaces prepared by 
fracture (involving crack propagation) , tend to be nearly 
self-affine fractal. Self-afiine fractal surfaces have mul- 
tiscale roughness, sometimes extending from the lateral 
size of the surface down to the atomic scale. A self-afhne 
fractal surface has the property that if part of the surface 
is magnified, with a magnification which in general is ap- 
propriately different in the direction perpendicular to the 
surface as compared to the lateral directions, the surface 
"looks the same" |2J| i.e., the statistical properties of the 
surface are invariant under this scale transformation. 

The most important property of a randomly rough sur- 
face is the surface roughness power spectrum defined as 



C{q) = 



1 



(27r) 



d'^x {h{x.)h{0))e''i-'- 



(3) 



Here /i(x) is surface height profile and (■ • ■) stands for 
ensemble average. We have assumed that the statisti- 
cal properties of the surface are translational invariant 
and isotropic so that C{q) depends only on the mag- 
nitude g = |q| of the wave- vector q. For a self-affine 
surface the power spectrum has the power-law behavior 
C{q) ^ q~'^^^^'^\ where the Hurst exponent H is re- 
lated to the fractal dimension D{ = 3 — H. Of course, 
for real surfaces this relation only holds in some finite 
wave- vector region qo < q < qi. Note that in many 
cases there is roll-off wave-vector qo below which C{q) 
is approximately constant. The mean of the square of 
the roughness profile can be obtained directly from C{q) 
using (T^ = (/i^W) = J d'^q C{q). 

For self-affine fractal surfaces r = A/Aq is uniquely 
determined by the root-mean-square (rms) roughness a 
and the fractal dimension Df. We have |2 71 



A/Ao 



dx {l + xi'^) 



l/2g-x 



(4) 



where f^ = j d^l (f'C{q). For the surfaces we use in our 
study, in Fig. ^ we show the ratio r = A/Aq both as 
a function of the root-mean-square roughness a, and as 
a function of Hurst exponent H. As expected, A/Aq 
increases with increasing rms-roughness, and decreas- 
ing Hurst exponent H (or increasing fractal dimension 
Df = 3 — H). Qualitatively, when Di increases at fixed 
rms-roughness, the short- wave length roughness increases 
while the long-wavelength roughness remains almost un- 
changed. 

We have used Molecular Dynamics calculations to 
study the influence of surface roughness on superhy- 
drophobicity. We have studied hydrocarbon liquid 
droplets on different self-affine fractal surfaces. The 
nano-droplet containes 2364 octane molecules CsHig at 
T = 300 K, which is between the melting and boiling 
points of octane. The fractal surfaces were generated as 
in Ref. |2g. Different fractal surfaces are obtained by 
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FIG. 1: The ratio A/Ao between the actual A and the nom- 
inal (or projected) Aq surface area, as a function of the root- 
mean-square roughness a when Hurst exponent H — 0.8, and 
as a function of Hurst exponent H for ct = 3 A. 



changing the root-mean-square roughness amplitude a, 
and the fractal dimension D[. The roll-off wave- vector 
for the rough surface is qo — 2t:/L with L = 38 A, and 
the magnitude of the short-distance cut-off wave-vector 
qi = tt/o, where a = 2.53 A is the substrate lattice con- 
stant. The (non-contact) cylindrical droplet diameter is 
about 104 A, and the size of the droplet-substrate con- 
tact area varies from w 115 A (case (a) in Fig. [2) to ~ 60 
A (case (c)). 

The lubricant molecules are described through the Op- 
timized Potential for Liquid Simulation (OPLS) 28, 2S|; 
this potential is known to provide density and viscosity 
of hydrocarbons close to the experimental one. We used 
the Lennard- Jones (L-J) interaction potential between 
droplet atoms and substrate atoms. The L-J parameters 
for a hydrophobic surface are chosen such that the Young 
contact angle is about 100° when a droplet sits on the flat 
surface. Because of the periodic boundary condition and 
the size of our system, the liquid droplet forms a cylinder 
with the central line along the y-axis, see Fig. |21 We fit 
the density profile of the droplet to a cylinder and obtain 
the contact angle 9 — 103° as indicated in Fig.|21for the 
droplet in contact with the flat substrate. 

The apparent contact angle, Oq, as a function of the 
root-mean-square roughness (rms), is shown in Fig. 01 
with the fractal dimension Df = 2.2. There is a strong 
increase in ^o with increasing rms-roughness amplitude. 
Fig. [S] shows how Og depends on the Hurst exponent 
H = 3 — D{. Note that 6*0 is almost independent of 
H. 

Accordingly to the Wenzel equation, the apparent con- 
tact angle 9o depends only on the surface roughness via 
the ratio r = A/Aq. Fig. ^ shows that as H decreases 
from 1 to 0.4 (i.e., D[ increases from 2 to 2.6), A/Aq 
increases by '^ 50%. However, the molecular dynamics 
calculations show that the apparent contact angle Oq is 
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FIG. 2: Snapshots for different root-mean-square roughness, 
(a) the droplet is in contact with the flat substrate, (b) and (c) 
are for rough substrates with the root-mean-square amplitude 
a = 2.3 A and a — 4.8 A, respectively. 
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FIG. 4: The contact angle as a function of the root-mean- 
square roughness a. The circle points are numerical results 
from the simulations, while the square points are obtained 
from the Cassie model (see Eq. I^J. Each data point is an 
average over several snap-shot configurations. The fractal di- 
mension is D{ = 2.2. 
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FIG. 5: The contact angle as a function of Hurst exponent 
H for the rms roughness ct = 3 A. The circles and squares have 
the same meaning as that in Fig. 0] The fractal dimension is 
Di = 3-H. 



FIG. 3: Determination of the contact angle 9 for the flat 
substrate. Side view. 



almost independent of the fractal dimension, see Fig. 
Thus the Wenzel equation cannot be used in the present 
situation. This is consistent with a visual inspection of 
the liquid-substrate interface which shows that on the 
rough substrates, the droplet is "riding" on the asperity 
top of the substrate, i.e., the droplet is in the Cassie state. 
In order to quantitatively verify this, we have calculated 
the distances h(x, y) between the bottom surface of the 
liquid drop and the rough substrate surface in the (ap- 
parent) contact area. From the distribution P{h) of these 
distances we obtain the fraction ip of the (projected) sur- 
face area where contact occurs: ip = L^ dhP{h), where 
hi is a cut-off distance to distinguish between contact 
and no-contact regions, which has to be comparable to 
the typical bond distance (we use hi = 4 A). Note that 



due to the thermal fluctuations ijj — "ipQ for flat surface is 
less than 1. Using the normalized (p = '0/V'Oi the Cassie 
model predicts the variation of the contact angle with a 
and H given in Fig. ^ and |5l (square points). 

Fig. ^ shows that the apparent contact angle 6*0 
increases strongly with increasing rms-roughness am- 
plitude, at fixed fractal dimension Df = 2.2, while 
it is nearly independent of the fractal dimension Df 
(see Fig. EJ. Since increasing the fractal dimension 
at constant rms roughness amplitude mainly increases 
the short-wavelength roughness, we conclude that the 
nanoscale wave length roughness doesn't matter so much 
in determining the contact angle for hydrophobic sur- 
faces, while the long wavelength roughness plays an im- 
portant role. We attribute this fact to the strong ther- 
mal fluctuations in the height (or width) h of the liquid- 
solid interface which occur on the nanoscale even for the 
flat substrate surface. Note also that in our model the 
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FIG. 6: The Advancing (circles) and receding (squares) con- 
tact angle 6^ as a function of time. The thermal equilibrium 
contact angle has been reached after ~ 1 ns, irrespectively of 
whether the initial contact angle is larger or smaller than the 
equilibrium angle. 



wall-wall interaction is long-ranged, decaying effectively 
as ^ l//i^, so there will be a contribution to the in- 
terfacial energy also for non-contacting surfaces which, 
of course, is not rigorously included in the macroscopic 
Cassie model. 

In Fig. we study the hysteresis [l^ in the contact an- 
gle 9. In one case a spherical droplet was attached to the 
substrate leading to a decrease in the contact angle with 
increasing time (advancing contact angle). In the other 
case the droplet was squeezed into a "pancake" -like shape 
by the upper wall and then released resulting in a con- 
tact angle which increases with time (receding contact 
angle). In both cases, the thermal equilibrium contact 
angle has been reached after ^ 1 ns. Thus, on macro- 
scopic time scales nano-scale roughness will not result in 
any hysteresis in the contact angle. This is in drastic 
contrast to simulation studies we have performed E30] for 
hydrophilic surfaces, where surface roughness results in 
strong pinning of the boundary line; for such surfaces it 
is therefore impossible to study (advancing or receding) 
droplet contact angles (as observed on macroscopic time 
scales) using molecular dynamics. 

Comparing the form of P{h) for the flat and the most 
rough surfaces shows that the system is in the Cassie 
state, but on the nanoscale the difference between the 
Cassie state and the Wenzel state is not so large due 
to the thermal fluctuations. This also explain why no 
hysteresis is observed: The Cassie state is the free energy 
minimum state and squeezing the droplet into a pancake 
shape does not push the system permanently into the 
Wenzel state because even if it would go into this state 
temporarily, the free energy barrier separating the Cassie 
and Wenzel states is so small that thermal fluctuations 
would almost instantaneously kick it back to the (free 
energy minimum) Cassie state. 

In most practical cases it is not possible to modify the 
surface roughness without simultaneously affecting the 



chemical nature of the surface. While this is obvious for 
crystalline materials, where surface roughening will re- 
sult in the exposure of new lattice planes with different 
intrinsic surface energy, it may also hold for amorphous- 
like materials, where surface roughening may result in 
a more open atomic surface structure, with an increased 
fraction of (weak) unsaturated bonds. In our model study 
a similar effect occurs, and some fraction of the change in 
the contact angle with increasing root-mean-square am- 
plitude may be associated with this effect. However, the 
most important result of our study, namely that the con- 
tact angle is mainly determined by the long-wavelength 
roughness, should not be affected by this fact. 

To summarize, we have studied the interaction between 
a liquid hydrocarbon nano-droplets and rough surfaces. 
The macroscopic contact angle 6*0 increases with increas- 
ing root-mean-square roughness amplitude a of the sur- 
face, but 00 is almost unchanged with increasing fractal 
dimension D{. There is almost no contact angle hystere- 
sis on the nanoscale. 
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